Growth activity during fingering in a porous Hele Shaw cell 
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We present in this paper an experimental study of the invasion activity during unstable drainage in 
a 2D random porous medium, when the (wetting) displaced fluid has a high viscosity with respect 
to that of the (non-wetting) displacing fluid, and for a range of almost two decades in capillary 
numbers corresponding to the transition between capillary and viscous fingering. We show that the 
invasion process takes place in an active zone within a characteristic screening length A from the 
tip of the most advanced finger. The invasion probability density is found to only depend on the 
distance z to the latter tip, and to be independent of the value for the capillary number C a - The 
mass density along the flow direction is related analytically to the invasion probability density, and 
the scaling with respect to the capillary number is consistent with a power law. Other quantities 
characteristic of the displacement process, such as the speed of the most advanced finger tip or the 
characteristic finger width, are also consistent with power laws of the capillary number. The link 
between the growth probability and the pressure field is studied analytically and an expression for 
the pressure in the defending fluid along the cluster is derived. The measured pressure are then 
compared with the corresponding simulated pressure field using this expression for the boundary 
condition on the cluster. 

PACS numbers: 47.20.Gv, 47.53.+n, 47.54.+r, 47.55.-t, 47.55.Mh, 68.05.-n, 68.05. Cf, 81.05.Rm. 
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I. INTRODUCTION 

Different types of unstable fluid displacements in 
porous media play an important role in many natural 
and commercial processes [lj, |2( . Developments of a bet- 
ter understanding of these processes therefore has a broad 
scientific interest as well as potentially huge economical 
benefits. The complex patterns observed in such pro- 
cesses have been extensively studied and modeled over 
the last decades, see 0, 0, y, E IH 111 an d references 
therein. 

The geometry of the displacement structures observed 
in immiscible two phase flow are in general controlled 
by the competition between viscous forces, gravitational 
forces, capillary forces; those various forces act on scales 
ranging from the pore scale to the system size. The rel- 
ative wettabilities, viscosities, and densities of the fluids, 
as well as the heterogeneity of the underlying porous me- 
dia, play an important role in the competition process. 
The relative magnitudes of viscous and capillary forces 
(on pore scale) arc quantified through the dimcnsionlcss 
capillary number C a — (fi w v / a 2 ) / (7 k) where fi w is the 
viscosity in the wetting (displaced) fluid, Vf is the fil- 
tration speed, a is the characteristic pore size, 7 is the 
interface tension, and k is the permeability of the porous 
medium. 

In this paper we address a drainage experiment in 
which non-wetting air displaces a high viscous wetting 
glycerin/water solution in a horizontal two-dimensional 
porous medium; hence, gravity has no influence on the 
displacement. We investigate the crossover regime be- 
tween the regime of slow displacement for which capillary 
forces control the dynamics of the invasion process and 



the geometry of the resulting invasion structure (capillary 
fingering), and that of fast displacements where viscous 
forces are dominant (viscous fingering) . We emphasize on 
the dependence of the invasion probability density cj>, or 
activity, on the distance to the most advanced finger tip 
along the interface. The invasion probability density <j) is 
the growth probability of the invasion structure; it is fun- 
damental because both the structure and the dynamics 
are controlled by this function. Growth probability has 
been discussed extensively in the past for DLA (Diffusion 
Limited Aggregation) simulations H H ES III 13 
where it was found to be a multifractal distribution or 
a harmonic measure 0, 0, 0, . A strong analogy 
exists between the structures obtained by DLA and by 
viscous fingering in a porous medium, as was first pointed 
out by Paterson 0] • However, the dynamics of the two 
processes differ in that there is no surface tension for the 
DLA, in contrast to drainage in porous media where sur- 
face tension gives rise to capillary pressure thresholds at 
the pore scale. The capillary pressure threshold values, 
introduce a lower cut off for the invasion probabilities, 
even for fast flows. In the slow displacement limit for 
which Vf ~ 0, the invasion process is entirely controlled 
by the fluctuations of the capillary threshold distribution 
inside the porous medium jy, Hq . 

Imbibition experiments (wetting fluid displacing a non 
wetting fluid) were previously performed in a quasi two- 
dimensional system 0], where the width of the 
viscous fingers was measured to scale with the capillary 
number as [riill8|| wj oc C~ 5 . This scaling relation was 
explained by a strong dynamic component of the capil- 
lary pressure We do not observe a strong dynamic 
component of the capillary pressure in our experiments 
(see below). The geometry of the invader for drainage is 
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also significantl y d ifferent from the invader structure of 
imbibition [3 El EES- 

In this study we study experimentally the growth prob- 
ability density <f>(z) as a function of the distance z (in 
the flow direction) from the most advanced finger tip, 
and its dependence on the extraction speed (or capillary 
number). We also investigate experimentally the mass 
density n(z) along the flow direction of the invader, and 
confront the behavior of the measured (f>(z) and n(z) to 
what we expect from analytical arguments. A calcula- 
tion of the z dependence of the pressure on the surface 
of the invader is presented, which yields the z depen- 
dence of the capillary pressure and shows a direct link 
with the measured growth probability density. Pressure 
measurements are performed in the model and compared 
with pressures simulated by solving the Laplace equation 
with this pressure boundary condition. Other features 
characteristic of the displacement, such as outermost tip 
velocity and the width of the invasion fingers, are also 
investigated. 

The present article is organized as follows. We first 
present the experimental method (section [nj . We then 
discuss the experimental results ( section ILTTjl . before con- 
cluding fsection HVI) . 



II. EXPERIMENTAL METHOD 

The experimental setup is shown in Fig. 2] The porous 
model consists of a mono-layer of glass beads of diameter 
a = 1 mm which is randomly spread between two contact 
papers pll |22T | . The model is a transparent rectangular 
box of dimensions L ■ W and thickness a. 

Two models of widths W — 430 mm and W = 215 mm 
have been used in the experiments; their other character- 
istics were identical. They are respectively referred to in 
the rest of the article as the "wide" and the "narrow" 
model. 

To prevent bending of the model a 2 cm thick glass 
plate and a 2 cm thick Plexiglas plate are placed on top 
of the model. To squeeze the beads and the contact pa- 
per together with the upper plate, a mylar membrane 
mounted on a 2.5 cm thick Plexiglas plate, below the 
model, is kept under a 3.5 m water pressure as a "pres- 
sure cushion" . The upper and the lower plates are kept 
together by clamps, and the side boundaries are sealed 
by a rectangular silicon rubber packing. Milled inlet and 
outlet channels are made in the upper Plexiglas plate. 
The distance between the inlet and outlet channels de- 
fine the length of the model L = 840 mm. One should 
also note that a few beads are removed from a small 
region near the center of the inlet channel, to initiate 
the invader in the center of the inlet. This is done to 
avoid edge effects appearing when the invader grows to 
the lateral boundaries of the model. The porosity of the 
models is measured to be 0.63 and the permeability is 
k = (0.0166 ± 0.0017) • 10~ 3 cm 2 = (1685 ± 175) Darcy. 

The defending wetting fluid used in all our experi- 
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FIG. 1: Sketch of the experimental setup with the light box for 
illumination, the porous model and the digital camera. The 
porous medium is sandwiched between two contact papers and 
kept together with a "pressure cushion". 



ments is a 90% — 10% by weight glycerol-water solution 
dyed with 0.1% Negrosine to increase the contrast be- 
tween the colored fluid and the invader. Air is used 
as the invading non-wetting fluid. The wetting glycerol- 
water solution has a viscosity of \x w ~ 0.165 Pa.s and 
a density of p w = 1235 kg.m~ 3 at room temperature. 
The corresponding parameters for the non-wetting air are 
Mnt« = 1.9T0~ 5 Pa.s and p nw — 116 kg. m -3 . The viscous 
ratio is thus M — fi nw /fi w ~ 10~ 4 . The surface tension 
between these two liquids is 7 = 6.4 • 10~ 2 N.m . The 
temperature in the defending fluid is controlled and mea- 
sured at the outlet of the model during each experiment, 
so as to accurately estimate the viscosity of the wetting 
fluid. 

The absolute pressure in the wetting liquid is mea- 
sured in the outlet channel and at a point at a distance 
of 280 mm (in the flow direction) from the inlet channel 
and 38 mm from the left boundary (looking in the flow di- 
rection) using Honeywell 26PCAFlow- Through pressure 
sensors. 

The invader is visualized by illuminating the model 
from below with a light box and taken pictures with a 
Kodak DCS 420 CCD camera, which is controlled by a 
computer over a SCSI connection. This computer records 
both the pictures and the pressure measurements. Each 
image contain 1536 x 1024 pixels, which corresponds to 
a spatial resolution of 0.55 mm per pixel or ~ 3.22 pixels 
per pore (1 mm 2 ); the color scale contains 256 Gray lev- 
els. The Gray level distribution of the image presents two 
peaks corresponding respectively to the white air-filled 
and dark gray glycerol-filled parts of the image. The im- 
age is filtered so as to obtain a clear boundary between 
the two phases, through a scheme that mainly consists 
in removing the background and thresholding at a gray 
level value between the two latter peaks. All further im- 
age treatments are performed on the resulting black and 
white image. 

To check possible dynamic components of the capillary 
pressure we performed gravity stabilized experiments by 
keeping the experimental model vertical |22| and extract- 
ing the glycerol/water mixture from the bottom of the 
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FIG. 2: Displacement structures obtained for different withdrawal rates: (a) C a = 0.027, (b) C a = 0.059 and (c) C a = 0.22. 
The images have been treated to separate the two phases, excluding the invasion structure close to the inlet. The black frame 
denotes the outer boundaries of the model, while the black spot close to the right edge of the model denotes the position of the 
pressure sensor. The simulated pressure field is shown superimposed on the image. Dark shadings correspond to low pressures 
while light shadings correspond to high pressures. 



model. The capillary pressure was measured by record- 
ing the pressure in the model as the stabilized fluid front 
approaches the sensor. No systematic dynamic effect on 
the capillary pressure was found. For the low injection 
rates the width of the fronts was further used to esti- 
mate the minimum and width W c of the capillary pres- 
sure threshold distribution. 

Throughout the paper the following coordinate system 
is used: (x, z') is the orthonormal frame describing the 
porous medium plane, with z' the spatial coordinate in 
the direction of the flow (positive in flow direction) . The 
position of the most advanced finger tip is denoted z' tip ; 
its speed along the z' axis is denoted v t i P — z' tip . The 
position along the z' axis computed with respect to that 
of the most advanced finger tip is z = z' tip — z' . Those 
coordinates are indicated in Fig. [21 and EJ 



III. RESULTS 

We present 12 experiments using the wide model for 
values of the capillary number C a ranging from 1.4 • 10~ 2 
to 3.6 • 10 _1 and 5 experiments using the narrow model, 
for capillary numbers ranging from 3 . 3 • 1 ~ 2 to 1.9-10. 
The latter series was conducted to check system size de- 
pendencies. For every experiment, we have carefully in- 
vestigated the invasion process. 



Fig-Gldisplays air clusters observed for the same porous 
medium, at three different flow rates. The "fingers" look 
visually thinner at increased invasion speed and more in- 
ternal trapping of the defender is observed at low speed. 
Hence, the displacement exhibits obvious capillary num- 
ber dependent features which will be discussed in details 
in part B below. In part A, we focus on the relation 
between the growth activity, the frozen structure left be- 
hind and the pressure field in front of the fingers. 



A. The relation between growth activity, frozen 
structure and fluid pressures. 

The growth activity has been investigated by measur- 
ing the growth probability density <f)(z) from series of 
images and performing pressure measurements. 



1. Definition of growth probability density 4>{z) and mass 
density n(z) . 

To investigate the growth process, images have been 
taken with constant time At between each image. The 
tip position of the longest finger is identified to find the 
coordinate system (x, z), and to be able to calculate the 
speed of the longest finger. The differential growth be- 
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FIG. 3: Two consecutive images taken at a time interval 
At ~ 15 during the experiment at C a = 0.059, drawn on 
top of each other. The invaded regions in the first image are 
painted light gray, the growth areas obtained by subtracting 
the first picture from the other one are painted black. The 
coordinate system used throughout the paper is also shown. 



tween two images is found by a direct image subtraction 
between two subsequent images. After the subtraction 
we typically have a collection of invaded pores represent- 
ing the growth (see Fig. |3J) . The growth density <j)*{z) 
is defined as the average number of filled pores within 
[z, z + Az] divided with Az. After an initial regime cor- 
responding to the time needed for the longest finger to 
propagate the order of the width of the porous medium 
W, 4>*{z) is found to be fairly independent of time up 
to a few percents variations. In a given experiment with 
constant C a , <j)*{z) is then averaged over all images ex- 
cluding this initial regime, to obtain a good average of 
the stationary growth function. The growth probability 
density 4>(z) = K <p*(z), where K is a normalization con- 
stant, is then found by normalizing 4>*{z) with respect to 
z such that 



(z) dz = 1 



(1) 



Note that in the reminder of the paper z is in units of 
pore size (a = 1 mm). 

The mass density of the frozen structure n{z) is defined 
as the average number of filled pores within [z, z + Az] 
divided with Az. The average is taken over all images in 
a given experiment with constant C Q . Both n(z) and 4>{z) 
appeared to be fairly robust with respect to the width Az 
of the analysis strips used to compute them. 
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FIG. 4: Plot of the position of the most advanced finger tip 
z' tip as a function of time. Data from the wide model. The 
time is rescaled by dividing with the break through time tb, 
which is the time the most advanced finger reach the outlet 
channel. The values of tb are 5347 s, 1204 s, 476 s and 256 s 
for the capillary numbers C' a equal to 0.027, 0.059, 0.12 and 
0.22 respectively. 



2. Growth activity and the frozen structure left behind 

For all experiments, the speed of the most advanced 
finger tip was observed to be fairly constant. Fig. 0] 
shows the position of the most advanced finger tip z' tip for 
different capillary numbers. After a short initial stage, 
the speed of the fingers saturates to a constant average 
value. Linear fits to the behavior z' tip as a function of 
time outside the initiation stage provide an average fin- 
ger tip speed vu p for all experiments. 

The measured invasion probability density function 
(j)(z) is plotted in Fig. [5] as a function of the distance to 
the finger tip for the two system sizes on a lin-log plot. 
An exponential like decay is seen for z/a < 100 with a 
deviation from exponential behavior for larger lengths. A 
characteristic decay length or "screening length" A is es- 
timated from linear fits to the lin-log data for z/a < 100 
(see Fig. E}. As we can see from these plots a nice 
data collapse is obtained, indicating that the invasion 
probability density <fi(z), and thus the screening length 
A, are independent of the capillary number for a given 
system. On the other hand, when comparing the two 
systems, the screening length A depends on the system 
size: A = (54 ± 10) mm for the wide model, and 
A = (34 ± 5) mm for the narrow model. The actual 
shape of function <p also seems to be weakly dependent 
on the system size. 

Relating the mass of the frozen structure, n(z), to the 
invasion probability density, 4>{z), and confronting the 
obtained relation to experimental results, provides new 
insights into the displacement process. The total number 
of invaded pores in a time interval [t, t + At] is R x At, 
where R is the number of invaded pores per time unit. 
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(a) z/a (b) z/a 

FIG. 5: Lin-log plot of the invasion probability density (j)(z) as a function of the distance z to the finger tip. The indicated 
distribution corresponds to (a) data from the wide experimental model and (b) from the narrow model. The solid lines in the 
curves corresponds to the model function ln(_Be _Z//A ) obtained from linear regression over \n((j>(z)) for z/a < 100. 



For a given flow rate Q, R is related to Q and to the 
characteristic pore volume V^, rc by the relation Q = Rx 
Vporc, so that 

W 

R=^v f , (2) 

where W is the width of the system, a is the characteristic 
pore size, and v/ is the Darcy or filtering velocity of the 
wetting fluid; for a given porous medium and fluid pair 
Vf oc C a . The number of invaded pores in the analysis 
strip defined by z £ [z, z + Az] during time interval [t, t + 
At] is then R At x <j>(z) Az. The tip position z' tip is further 
given by z' tip (t) = z' Q + Vf lp t where «tip is the speed of the 
finger tip (assumed to be constant) and z' = z' tip (t = 0). 
This is a fairly good approximation after a short initial 
regime as seen in Fig. 

The total number of invaded pores in an analysis strip 
at distance z from the finger tip and in a strip of width 
Az is thus given by: 

n(z) Az = R ( 4>{z{t')] Az dt' , (3) 

J to 

where to is the time at which z' tip = z' . Taking advantage 
of the linear relation between the coordinate z and tip 
speed wtip, z(t) = z' tip (t) - z[ ip (t ) = (t-t ) v tip , Eq. © 
becomes 

n(z) = — [ <t>(z) dz = — $(z) . (4) 

"tip Jo v tip 

Using Eq. (J2J) we finally obtain the relation 

n(z) = $(z) oc — $(z) , (5) 

a va v «tip 

which relates the linear density of invaded pores (or 
"cluster-mass" density), n{z), to the cumulative invasion 
probability density distribution, <&(z). This relation is 



confirmed by Fig. El where n(z) x (a 2 v t i P / (Wvf j) is plot- 
ted as a function of z/a. All experimental plots collapse, 
confirming that there is one single cumulative probability 
distribution Q(z) for the system for all experiments at dif- 
ferent extraction speeds. The function Q(z), computed 
as an average function from all cumulative probability 
functions for the various experiments, is plotted in Fig. El 
as a plain line. 

The insert of Fig. ©shows {a 2 v t i V )/{Wvf){n ocl — n{z)) 
on a lin-log scale, where rioo = Wvf /a 2 v t i p . The solid 
lines represents (1 — $(z)) and the dashed lines (1 — e~ z / A ) 
which would be the model function for a pure exponential 
4> (A is the screening length evaluated before). 

From the results presented above we conclude that the 
active invasion zone is defined by a screening length A 
which is constant for a given porous media and liquid 
pair and at a range of capillary numbers of two decades. 
However, we expect this result to be valid only for suffi- 
ciently high filtration speed Vf. Indeed, on the one hand, 
the capillary fingering regime (C a ~ 0) corresponds to 
an invasion that is controlled by fluctuations in the cap- 
illary threshold pressures, so that invasion occurs along 
the whole front There is no well defined finger tip 

or growth direction in that limit. The width of the capil- 
lary threshold pressure distribution W c is larger than the 
viscous pressure drop over the whole system and defining 
a screening length or active zone is not meaningful. 

When the length of the system is larger than its width, 
it is found from both pressure measurements and simu- 
lations that the decay in pressure into the structure from 
the longest finger occurs on a length scale of the order of 
the width of the system (see Fig.|5J). We therefore expect 
W and not L to be the relevant length scale for the decay 
of the pressure field close to the tip. Viscous forces can 
therefore be considered to dominate capillary pressures 
if the following criteria is met: 

W C <^1 , (6) 

K 
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(a) z/a (b) z/a 

FIG. 6: Plot of the rescaled average mass density of non wetting fluid A x n(z) where A — a 2 v t i P /(W Vf) inside the model as 
a function of the distance to the tip of the most advanced finger z: (a) the wide model and (b) the narrow model. The average 
cumulative invasion probability function <&(z) is plotted on top of the curves for comparison in the two cases. The inserts of 
(a) and (b) show the plots In [l — {a? vup/(W «/)) n(z)] = In [(a 2 vu p )/(W + In [{noo — n { z ))] as a function of z. The solid 
lines in the inserts correspond to ln(l — $(z)), and the dashed lines has the slope 1/A where A is the screening length of 4>{z) 
found from Fig. [3| i.e. respectively A = (54 ± 10) mm and (34 ± 5) for the wide model and narrow model. 



or if we assume W c ~ P c (which is the case here) 

C a >a/W . (7) 

For our system a/W ~ 10 -3 , which is an order of mag- 
nitude smaller than our lowest capillary number. 

On the other hand, for situations where the "pure 
viscous fingering" in a random porous media has been 
reached, there is no trapping of wetting liquid inside the 
fingers, which reached the lower one pore width limit (at 
C a « 0.2 in our system). Whether the screening length 
or active zone have the same width or behave identically 
as for lower capillary numbers is not clear. We believe 
that the screening by the most advanced finger is a vis- 
cous effect, which remains important as the displacement 
speed increases. In this one pore limit however, the tip 
speed dependence on the capillary number is modified, 
as will be further detailed in section ITlI Bl 

3. The relation between the growth probability density <f>{z) 
and the fluid pressure. 

Fig. E| shows the dependence of the pressure difference 
AP(z) — P{z) — -P(oo) in the wetting liquid as a func- 
tion of the distance z to the outermost tip for different 
capillary numbers. Here P(oo) = Pq — P c (oo) is the liq- 
uid pressure along the interface far behind the finger tip, 
with Pq the pressure in the non- wetting liquid and P c (oo) 
the capillary pressure in this stagnant zone. It is impor- 
tant to note that the pressure P(z) is measured on the 
side of the model (indicated in Fig. [2J while the fingers 
are propagating in the central part of the model (Fig.|2J). 

The pressure seems to be linearly dependent on the 
distance from the tip during a first stage before the tip 



reaches the sensor. In a second stage, after the finger tip 
has passed the sensor, pressure relaxes and reaches the 
value P(oo). A closer inspection of the pressure curves 
(see Fig. [7J shows that there is no clear systematic de- 
pendence of the pressure relaxation on the capillary num- 
ber. The pressure difference P(z) — P(oo) decays with 
approximately the same length for the different capillary 
numbers (see insert of FigEJ. This indicates that the de- 
tails of the internal structure of the "fingers" do not have 
a strong influence on the pressure field on large scales. 

The pressure measurements are related to the invasion 
activity by the following considerations. Let us consider 
the local speed of an interface located in an arbitrary 
pore throat between two pores, one filled with air and the 
other with the wetting-liquid. Let P(x, z) be the pressure 
in the wetting liquid and Pt(x,z) be the capillary pres- 
sure threshold value to invade that pore. Note that this 
is different from the pressure P(z) defined as the pressure 
measured on the side of the model at the sensor position. 
The pore throat at position (x, z) is passed under the con- 
dition that the capillary pressure P c (x, z) = Pq — P(x, z) 
is larger than the capillary threshold pressure Pt(x, z) at 
this position. If invasion occurs, a characteristic value of 
the speed of the interface will be: 

fi a 

Let N(Pt(x, z)) be the capillary pressure distribution. 
For the sake of simplicity, we assume a fiat capillary pres- 
sure distribution with lower limit P t mm , upper limit P t max 
and width W c . Then, the expectational value of the inter- 
face velocity (average value over the capillary threshold 
distribution), while the pore is getting invaded will be 



7 



• >»o 



-0.5 



a -i 

Oh 

< 



-1.5 



T 



0*4*" 



- w 



(a) 



-100 



1 - 3 
a. 

S-4 
c 

- -5 
-6 





I 

100 





c 

a 


= 0.027 


V 


c 

a 


= 0.059 


+ 


c 

a 


= 0.12 


* 


c 

a 


= 0.22 



I ()() 200 300 400 

j I i I l_ 

200 300 
z/a 



500 

I 

400 



500 



-0.25 



-0.5 
£-0.75 

Oh 

-1 



_6> 

50 



-1.25 
-1.5 



< 



o 


C 

a 


= 0.0334 


o 


c 

a 


= 0.079 


< 


c = 


= 0.19 



-100 



00 




400 



500 



FIG. 7: PZoi o/ i/ie pressure difference AP(z) — P(z) — P(oo) measured by the pressure sensor located in the position (x 3 , z' s ) = 
(38 mm, 280 mm). The insert is the same data plotted in a semi-log plot where AP(z) is scaled with AP(0) = P(z — 0) — P (oo) 
to illustrate the decay of the pressure field. Figure (a) is for the wide model and Figure (b) is for the narrow model. 



(v(x,z)) = 



1 



Po - P(x, z) 



P -P(x,z) 2 

— (P - P(x, z) - P t (x, z)) dP t = —{P - P(x, z) - P t min ) . (9) 



(if i 



(If l: 



r 



Here, P t mm is the minimum of the distribution for capil- 
lary threshold; when Po ~ P{x, z) goes to that minimum, 
the expectational value for the speed of the interface goes 
to zero. The growth probability density 4>{x, z) for the 
invasion structure within a time [t, t + At] at a position 
(x, z) is proportional to {v{x, z)) times the probability 
p(x, z) that the throat gets invaded; hence 



(x,z) = C (v(x,z))p(x,z) 



(10) 



where C is a normalization constant, which we can find 
by integrating the above equation along the invasion front 
S: 



Is 



(x,z)dl — C / (v(x, z) ) p(x, z) dl 
Js 

C f C 
1 = — / a (v(x, z) } p(x, z)dl — — Q 
a Js a 



where Q is the flow rate, thus: 
C-- 

C Q 



(ii) 

(12) 
(13) 



Since we have assumed a flat capillary threshold dis- 
tribution of width W c , the probability that the pore at 
position {x, z) gets invaded is 



P(x,y) = ±r (P -P(X,Z)-P? 



(14) 



From Eq. JH1 - I14|l we obtain for the growth probability 
density 4>{x, z) in position {x, z): 



(x,z) 



QfiW c 



r>min 



P(x,z)Y 



(15) 



Averaging this expression over x and introducing the 
number of interface sites at a distance z from the tip, 
/(z), we obtain the invasion probability density 4>(z) as 



(*) = /(*) {4>{x,z)) x 



(16) 



(z) = f(z) 



QfiW c 



(p Q -pr n -{p(x lZ )) x ) 2 , a?) 



for which we have assumed that P(x, z) is a function of 
z only (lowest order approximation). Eq. (|17[) yields 



(P(x,z))=P 



pnun 



cj>{z) QW c fi 

/(*) * 



(18) 



which can be rewritten by introducing the relation be- 
tween the flow rate and the capillary number. Accord- 
ingly the average pressure in the wetting fluid in the im- 
mediate vicinity of the interface and at position z is re- 
lated to the activity 4>{z) according to 



(P(x,z)) x =P - 



pmin 



c al w c 



w 



a f(z) 



(19) 



Let us now look closer at the "snapshots" of the ex- 
periments shown in Fig. [3 For z — 0, the last correc- 
tion term in Eq. 1)19(1 is 170 Pa for the fastest experi- 
ment (C a = 0.22), and 65 Pa for the slowest experiments 
(C a = 0.027). At the same moment, the imposed ex- 
ternal pressures in the outlet channel are 3055 Pa for 
the fastest and 625 Pa for the slowest experiments. The 
minimum capillary pressure is estimated to 373 Pa, and 
the width of the capillary distribution to 200 Pa. This 
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indicates that the correction term in Eq. i|I9[l should 
not be neglected. In Fig. [21 is shown the gray scale 
map of the pressure field at a particular time, simu- 
lated from the displacement structures obtained exper- 
imentally. A very strong screening is seen for all in- 
jection rates. The large scale structure of the pressure 
field in the vicinity of tip of the longest finger looks visu- 
ally very similar even if the invader structure is quite 
different. In the simulations of the pressure field we 
have used Eq. lj*H?|l to set the proper boundary condi- 
tions. The pressure field has been calculated by solving 
the Laplace equation for the pressure using a conjugate 
gradient method [23J. We used the boundary condition 
given by Eq. (|19|l on the cluster and the inlet line. As 
boundary condition on the outlet we used the pressure 
P(oo) — AP tot , where — AP to t is the total viscous pres- 
sure drop imposed in the corresponding experiment at 
that moment. To obtain P(oo) = Pq — P c (oo), the cap- 
illary pressure P c (oo) was measured in the experiments 
for large values of z. Fig. [S] shows the simulated pres- 
sure AP(x s ,z) = P(x s ,z) — P(oo) as a function of the 
z-coordinate relative to the tip position defined as previ- 
ously, at a fixed lateral position x s corresponding to the 
^-coordinate of the pressure sensor. It is important to 
note that this is somewhat different from the experiments 
since the pressure is measured at different z' positions, 
but at the same time, i.e. with a fixed geometry of the 
invasion cluster, while in the experiments the pressure is 
measured at a fixed z' position, at different times corre- 
sponding to various stages of the invasion cluster. The 
length scale of the decay of the pressure for z > is very 
similar in the experiments and the simulations (see com- 
ment below). However the pressure difference AP(z) in 
the simulations is lower than the AP(z) measured by the 
sensor (at position (x s ,z' s ) = (38 mm, 280 mm)) in the 
model (see Fig. This is due to the strong boundary 
effects of the pressure close to the outlet channel: the tips 
in the simulations situations are very close to this bound- 
ary along which the pressure is fixed (see Fig. |2J), while 
in the situations corresponding to the measurements at 
z < plotted in Fig.[7| the outlet boundary was far ahead 
of the finger tip, and the pressure boundary condition was 
equivalent to an imposed gradient at infinite distance. To 
check the importance of this boundary effect on the mag- 
nitude of the pressure difference, the simulated pressure 
has been compared with the pressure difference AP(z) 
evaluated from measurements at the outlet channel, as 
the finger tip progressed further than the stage corre- 
sponding to the simulations. The agreement in Fig. [S] 
between the simulated pressure and the data points cor- 
responding to the outlet channel measurements is then 
satisfactory. 

To compare the length scale of the decay of the pres- 
sure for z > between the simulations and the experi- 
ment, we then compare the pressure data measured in- 
side the model at sensor position, to the simulation data 
scaled by a factor such as AP(x s , z = 0) would be equal 
in experiments and simulations. Such rescaling of the 




500 



FIG. 8: The simulated pressure along the line with x s = 
38 mm ( corresponding to the x coordinate of the pressure sen- 
sor) for the invasion structures in Fig. The data points in 
the main graph show the corresponding pressures measured at 
the outlet sensor. In the insert we plot scaled simulation data 
with corresponding pressure data measured inside the model 
at (x s ,z' s ) — (38 mm, 280 mm). 



simulation pressure profile simply corresponds to the re- 
sult of an identical simulation still carried on the invasion 
clusters of Fig. with identical boundary conditions de- 
rived from the growth density function for the pressure 
along the clusters, but where the imposed pressure along 
the bottom boundary is such that the pressure at point 
{x s ,z' tip ) would coincide with the pressure measured in 
the experiments when the tip passed at the same height 
as the sensor, i.e. when z' tip = z' s . This ensures that the 
pressure gradient and pressure value in the region around 
the tip of the invading cluster are of the same order in 
these rescaled simulations and in the experimental stages 
corresponding to z ~ in Fig. |7| which is a first order 
technique to correct for the strong boundary effect and 
compare with these experimental situations where the 
bottom boundary is much further away. The pressure 
measured in the experiments at sensor position and this 
scaled simulation data are plotted in the insert of Fig. |S] 
This comparison shows that the decay in the pressure 
happens at comparable length scales in the simulations 
and experiments. 

Eventually, the local structure of the finger and the 
lateral x-distance from the invader to the pressure sen- 
sor will also have an important influence on the pressure 
field. The difference in the pressure field between the left 
and the right side of the finger (looking in the flow direc- 
tion) in Fig. [2] illustrates this point. The deviation be- 
tween the experimental data points and the simulations 
for the lowest capillary number of the main part of Fig. [S] 
may be explained by this effect. As the lateral position 
xup of the invading structure moves during the experi- 
ment, and is importantly varying from an experiment to 
the next, this effect also explains the important disper- 
sion of the scaled pressure drops AP(z)/AP(0) observed 
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in the insert of Fig. (b) . 



B. Capillary number dependent features 

As stated in the introduction to section II I II Fig. [21 
clearly shows that some features of the invading clus- 
ter depend on the capillary number. The mass n(z) of 
the invasion cluster obviously decreases with increasing 
capillary number; in relation to this, the speed of the 
most advanced finger tip, Vf lp , increases with the capil- 
lary number, and there is a systematic trend for fingers 
to become thinner as capillary number increases. In the 
following we first present results relative to the "mass 
density" in the stagnant zone, n^, and to the velocity of 
the most advanced fingertip, Utip- in the end we discuss 
the results relative to measurements of the characteristic 
width of the finger- like structures, the definition of which 
are not as straightforward and clear as those of rioo and 

Vtip- 

The evolution of the average mass density in the stag- 
nant zone Hoc as a function of the capillary number is 
presented in Fig. 0(a) on a log- log scale. The data is con- 
sistent with a scaling law in the form cx C~ Q , with a 
scaling exponent a — 0.65 ± 0.05 for both the wide and 
narrow models. Here has been measured by fitting 
the function [1 — exp(— z/X)] with both parameters 
free to our measured n(z) data. Due to that dependence 
of the mass of the invasion cluster on the capillary num- 
ber (Fig. (a)), the speed of the most advanced finger 
tip, vup, is expected to depend on the filtration speed 
or capillary number in a non-linear way. The saturated 
mass density and the speed of the most advanced finger 
are related to each other through Eq. |SJ according to 

C 

vtip oc — . (20) 

77-00 

Based on that argument, C a /vti p should therefore scale in 
the same way as with respect to the capillary number. 
In Fig.El(b), the quantity C a /vti p is plotted as a function 
of C a using a log- log scale. The plot is consistent with the 
expected scaling (12011 and the result for the mass density 
presented above. 

The study of the dependency of the finger width on 
the capillary number is somewhat less straightforward, 
because our invading clusters structures exhibit exten- 
sive branching and display "fingers" both at small scales 
as "capillary fingers" and at large scales as "viscous fin- 
gers" . Thus, a precise definition of a finger, and further- 
more a finger width, is not an easy task for those struc- 
tures. A possible method to determine the viscous finger 
width would consist in finding the characteristic crossover 
length between geometric features characteristic of vis- 
cous fingering and those characteristic of capillary finger- 
ing from the density-density correlation function of the 
structures. However, due to the small difference in fractal 
dimension between the two regimes, 1.83 ± 0.01 || 0] 



for capillary fingering and 1.62 ± 0.04 [2lt |24| for vis- 
cous fingering, larger systems would be necessary for this 
method to be accurate enough. An experimental deter- 
mination of the characteristic width Wf for viscous fin- 
gers was previously obtained for imbibition experiments 
|17|. for which the characteristic finger width can be de- 
fined and found in a more straightforward manner. The 
obtained scaling was Wf cx C~ ' 5 . In those imbibition 
experiments, the finger width Wf was measured as the 
average length of cut-segments perpendicular to the flow 
direction. This method can also be applied in our ex- 
periments, but due to the small scale fractal nature of 
the invasion front, trapping of wetting fluid inside the 
fingers and capillary fingering on small length scales, it 
is not obvious which length scales are being probed with 
this method. The results that we obtain are plotted as 
a function of the capillary number in Fig. 1101 Clusters 
of wetting liquid trapped behind the displacement front 
have been removed from the picture prior to analysis. 
We then define the front width wj as the average over 
z and time of the length of the intersects between the 
invasion cluster emptied from these trapped regions and 
cuts perpendicular to the flow direction. 

The measurements are consistent with a scaling law in 
the form w f cx C' 13 , with (3 = 0.75 ± 0.05. This is signif- 
icantly different from what was measured for imbibition. 
It also differs significantly from the scaling law expected 
from theoretical arguments for percolation in a desta- 
bilizing gradient |25l |26| for two-dimensional systems: 
Wf oc C-P with = 0.57. In our experiments, the desta- 
bilizing field (pressure) is highly inhomogeneous, which 
may explains why the behavior expected from the perco- 
lation in a gradient theory is not really observed. 

From Fig. and ^3 the observed scalings appear to 
be valid for a limited range of capillary numbers. For 
high capillary numbers the observed scaling breaks down 
for C a ~ 0.2, which corresponds to situations where the 
characteristic finger width have reached the one pore 
limit. At the other limit, for small capillary numbers 
we don't seem to reach the lower limit in capillary num- 
ber. But we expect that the observed scaling brake down 
for capillary numbers smaller than the criteria given in 
Eq. O, C a ~ 10- 3 . 



IV. CONCLUSION AND PROSPECTS 

We have studied the dynamics of the invasion process 
observed during drainage in a two-dimensional porous 
medium, for extraction speeds that result in an unsta- 
ble fingering of the displacing non- wetting fluid into the 
displaced wetting fluid. 

Our main finding is that for a given porous medium, 
the displacement is controlled by an invasion probability 
density that only depends on the distance of the point 
where it is measured to the tip of the most advanced 
finger tip, and is independent of the capillary number. 
The decay of this invasion probability density, 4>(z), de- 
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FIG. 9: Log-log plots of (a) the saturated mass density, rioo, and (b) of the speed of the most advanced finger , vu p , as a function 
of the capillary number, for the two set of experiments. Both plots are consistent with a scaling in the form oc C a /vu P oc C~ a , 
with a = 0.65 ± 0.05. 



o wide model 




log 10 (C) 

FIG. 10: Log-log plot of the measured characteristic width of 
the finger like structures as a function of the capillary number, 
for both the wide and narrow model. The data is consistent 
with a scaling of the finger width in the form u>f oc Ca 13 , with 
(5 = 0.75 ± 0.05. 

fines an active zone for the invasion process, outside of 
which the viscous pressure field can be considered to be 
screened by the invasion structure. In particular, parts 
of the invasion structure lying outside this active zone 
are frozen and do not evolve in time any more. The 
size of the active zone, of characteristic screening length, 
A, was found to be independent of the capillary number 
for a wide range of injection rates. In addition, exper- 



iments carried out on models with two different widths 
suggested that the invasion probability density appears 
to be capillary number independent, its actual shape be- 
ing possibly fixed by the system size. While the invasion 
process is described by an invasion probability density 
that is independent of the capillary number, the invasion 
speed and displaced volume in the stagnant zone were 
found to scale on the capillary according to power laws, 
"oc oc C a /v tip oc C-°' 65 . 

The link between the growth probability and the pres- 
sure field has been studied. An expression for the pres- 
sure boundary condition on the cluster has been calcu- 
lated which relates the pressure on the interface of the 
invader to the growth probability density function <f>(z). 
The measured pressure has been compared to the corre- 
sponding simulated pressure by solving the Laplace equa- 
tion for the pressure field using this expression for the 
boundary condition on the cluster. A good agreement is 
found between the simulations and the experiments. 

System size dependencies should be subject to further 
investigations, both experimentally and by means of com- 
puter simulations. 
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